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ABSTRACT (Continued) 


Several simulations are used to evaluate the algorithm. Estimated natural 
frequencies and damping ratios are compared with simulation values. The 
algorithm is also applied to wind tunnel data in an off-line mode. The 
results are used to develop preliminary guidelines for effective use of 
the algorithm. 
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CHAPTER I 

INTRODUCTION AND SUMMARY 


1.1 INTRODUCTION 

The XV- IS is an advanced tilt rotor research aircraft 
(Figure 1-1) currently under development by the National Aero- 
nautics and Space Administration and the U.S. Army Air Mobility 
R§D Laboratory. An important element of the aircraft develop- 
ment is an extensive full-scale wind tunnel test to be conducted 
in the Ames Research Center 40- by 80-foot wind tunnel. Since 
wind tunnel scale model and full-scale tests of the tilt rotor 
concept have been effective in previous years, it is desired 
to conduct the forthcoming XV-15 tests in a comprehensive 
manner and to continue to minimize uncertainties in system 
characteristics . A previous study evaluated specific test 
requirements which impact the conduct of these tests [1] . 



Figure 1.1 The NASA/Army XV-15 
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The present study describes a parameter identification 
method by which the dynamic stability characteristics of the 
XV-15 may be calculated on-line, during a wind tunnel test run. 
This is accomplished through the use of an advanced parameter 
identification algorithm, developed specifically for this 
application. The development of this algorithm and its appli- 
cation to several simulations and cantilever wing wind tunnel 
test data are the subjects of this report. 


1 . 2 BACKGROUND 


It is now recognized that the use of comprehensive soft- 
ware for aircraft wind tunnel and flight test data reduction 
can significantly enhance vehicle usefulness by producing 
the following results: 

(1) improvement in knowledge of vehicle static and dynamic 

characteristics : comprehensive software, accounting 

for typical disturbances in the data, can provide 
more accurate definition of the desired test objec- 
tive parameters (e.g. stability); 

(2) reduction of test costs : more engineering results 

can be obtained with less test and analysis time and 
specialized personnel; and 

(3) enhancement of test safety : for real time software, 

rapid determination of test system status allows 
faster response to problems. 

The testing of advanced rotorcraft, such as the XV-15, is 
a current requirement for such software. This vehicle sets a 
particularly stringent data reduction requirement for many 
reasons including the following: 

(1) Previous tunnel experience of a similar configuration , 
the XV- 3, demonstrated fundamental phenomena whose 
subsequent successful analysis has led to the cur- 
rent XV-15 design [2-8]. This experience indicates 
the importance of the analysis of the tunnel test 
data analysis. 
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( 2 ) 


Extensive tunnel testing of the XV- 15 is planned to 
thoroughly investigate regimes of the subsequent 
flight tests, and also provide a basis for evaluation 
of the flight test data reduction system. This 
requires a data reduction capability which is not 
only comprehensive (e.g. applicable to many flight 
phenomena) , but may be applicable to flight data 
as well. 

(3) The advanced configuration of the XV- 15 requires an 
advanced multivariable data reduction technique for 
on-line estimation of modal frequency and damping. 

This is a requirement for a new technique, compatible 
with existing Ames 40- by 80-foot wind tunnel 
computers . 

Several previous analytical developments have been made 
to aid the wind tunnel tests. Some of the important develop- 
ments are: 

(1) a comprehensive XV- 15 analysis has been developed and 
used for study of partial and full-scale models [5-8], 
simulation of transient responses [1], and linear 
stability analysis [1,5-8]; 

(2) investigations of input and instrumentation require - 

ments for the XV- 15 tunnel test, based on the model 
of (1) [1] ; and 

(3) study of various off-line data reduction techniques 

for determining stability estimates from tunnel test 
data, including: (a) time series analysis techniques, 

and (b) maximum likelihood identification techniques 

[ 1 ]. 

To complete the XV-15 tunnel test requirement, it was 
considered necessary to develop a new algorithm to extract 
the most information from the wind tunnel test data. 


1.3 METHOD OF APPROACH 

The requirements of the previous section demand a data 
processing algorithm with many special features. The algorithm 
should produce accurate estimates of damping ratios and fre- 
quencies of various lightly damped modes, should require small 
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computation time and storage, should account for measure- 
ment and process noise, and should require minimum a priori 
information. 

To meet these requirements, the dynamic model of the 
XV-15 tilt rotor is written in a continuous time, autoregres- 
sive moving - average (ARMA) form. The ARMA representation is 
equivalent to a state variable form or a generalized state 
form (see Appendix A) . 

The estimates of the ARMA model parameters are obtained 
using a time domain instrumental variable technique. The 
instrumental variables technique minimizes the effects of 
process and measurement noise and provides a computationally 
efficient technique which does not require any a priori esti- 
mates. This time domain formulation is converted into a fre- 
quency domain formulation for ease of implementation. 


1.4 PRINCIPAL CONTRIBUTIONS 


The contributions of this work include basic theoretical 
developments and applications study. The theoretical develop- 
ments lead to a modal estimation algorithm and the application 
results develop procedures for using the algorithm effectively. 

New theoretical contributions include the following: 

(1) A special autoregressive moving average (ARMA) form 
is developed for dynamic systems, with general model 
structure, for the purpose of measuring dynamic 
stability. The coefficients of the system are 
directly obtained from the autoregressive part of 
the model. 

(2) An instrumental variables algorithm is developed 
for estimating the parameters of the ARMA model. 

The identification procedure gives the transfer 
function numerator and denominator polynomials 
directly. The algorithm can use a multi- input, 
multi-output data efficiently and minimizes the 
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effects of measurement and process noise. The fre- 
quency range over which the model is required may 
be specified. 

( 3 ) A heuristic approach for model order determination 
is developed based on approximate cancellation of 
identified poles and zeros. 

The application study has lead to the following major 
contributions: 

(1) Tentative procedures have been developed for set- 
ting various algorithm parameters for the most 
effective use of the algorithm. 

(2) Procedures for obtaining good numerical conditioning 
in the application of the algorithm, particularly 
with many degrees -of- freedom. 

These contributions have established the superiority of 
this system identification-based algorithm over classical 
approaches. Further work is required, however, in maximizing 
its effectiveness in on-line application to measure stability 
in the Ames 40- by 80-foot Wind Tunnel. 

1 . 5 SUMMARY 

Subsequent chapters of this report are organized as 
follows. Chapter II describes qualitatively approaches to 
estimating stability parameters of aeroelastic modes. Chapter 
III discusses the problem, the choice of models, and the 
methodology used in the development of the parameter estimation 
algorithm. The algorithm is applied to ground resonance and 
cantilever wing simulations . The algorithm is further studied 
using XV-15 simulations in Chapter IV. Chapter V describes 
the application of the algorithm to a scaled-model wind tunnel 
test. The conclusions are given in Chapter VI. 
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Appendix A describes the relationship between various 
forms of dynamic models. A detailed derivation of the algo- 
rithm is given in Appendix B. Appendix C shows the salient 
features of the XV-15 simulation used in Chapter IV. 
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CHAPTER II 

APPROACHES TO DYNAMIC STABILITY MEASUREMENT 


2 . 1 INTRODUCTION 

Interest in damping ratios and natural frequencies of 
modes is as old as linear dynamic models. Damping ratios 
determine the extent of stability and the natural frequencies 
determine the "speed” of a system. Early work in the deter- 
mination of dynamic stability parameters from test data xvas 
performed to evaluate handling qualities of airplanes. This 
work can be traced back to Milliken [32], who made important 
observations not only from the viewpoint of data processing 
but test design as well. Prony's response function methods, 
and -analog matching techniques, were used extensively with 
limited success . The development of high speed digital com- 
puters in the early 1960's considerably advanced the computa- 
tional capabilities without a parallel development in estima- 
tion techniques. Only in the last five years has sufficient 
work been done to enable accurate estimation of dynamic stabil 
ity parameters from limited noisy test data. A qualitative 
review of current approaches to the problem of determining 
stability (i.e., modal damping ratio) from experimental meas- 
urements is presented in the following, as background to the 
new parameter identification algorithm discussed in the 
remainder of this report. 

2.2 PREVIOUS APPROACHES 

Estimates of modal frequency and damping ratio have been 
obtained using well -understood methods of spectral analysis 
and known dynamic characteristics of second-order systems 
[10-13] . These methods generally give good results: O) when 
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modes are not closely coupled, (b) when noise levels are low 
(i.e., when the second-order characteristics of a mode are 
discernible), and (c) with subjective interpretation of an 
experienced analyst. Several methods that have given accept- 
able results are discussed in Ref. 9 and are reviewed in the 
following. These methods include: (1) transient decay, 

(2) moving block analysis of transient decay, (3) the method 
of random decrement signatures; (4) spectral analysis, and 
(5) transfer function analysis. 

Transient decay methods [13] determine damping ratio from 
the envelope of a decaying oscillation. If the oscillation 
is assumed to be second-order, damping ratio may be deter- 
mined from response amplitudes at two points, one or more 
cycles apart. This method is susceptible to error from noise 
and extraneous modes which cause the appearance of the signal 
to deviate from that of a second-order system. In some cases, 
filters can be used to make the mode of interest more promin- 
ent, and careful input selection can often avoid exciting 
other modes . 

In the method of "moving block" analysis, portions of 
the response signal are Fourier-analyzed at successive points 
in time. The damping ratio is then calculated as in the 
transient decay method, but using the magnitude of the spectral 
line at the fundamental frequency of the mode of interest 
rather than the measured signal amplitude. This method has 
the advantage of filtering out a portion of noise, if there 
is a reasonable frequency separation between the noise and 
the signal. Damping ratio calculation may be mechanized rela- 
tively easily, an advantage over the transient decay method. 

In practice, a fast Fourier transform (FFT) method may be 
used, but three or four periods should be included in each 
block to achieve sufficient accuracy. The method is thus 
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applicable only to low damped modes such as are encountered 
in aeroelastics . 

The method of random decrement (RANDOMDEC) signatures 
[17,18] is a technique for determining the impulse response 
of a linear system from measurements of its response to random 
disturbances. When the system is under random excitation, 
as from wind tunnel turbulence or a specific random input, 
measurements are taken at times which have identical initial 
conditions on the system states. In taking repeated measure- 
ments in this manner, random effects cancel out, leaving in 
principle the system response to initial conditions. The 
method as developed by Cole [17] may be easily mechanized on 
a digital computer. The drawbacks of the method are: 

(a) the presence of coupling modes reduces the effectiveness 
of the technique, (b) it is difficult to achieve the same 
initial condition on every state variable (not only on 
every output] in each run, (c) most of the data cannot be 
used because of the requirement of the same initial condi- 
tion, therefore the test time is too high, and (d) the 
errors in triggering are not averaged out. 

Analysis of the power spectrum of the measured signal 
can be used to estimate the natural frequency of a mode by 
the location of its resonant peak. For this purpose, the sys- 
tem exciting function spectrum should be relatively "flat” 
near the resonant frequency [14]. By relating measured 
spectra, which may be averaged over successive measurements, 
to the spectrum of a second-order system, estimates of damping 
ratio can be obtained. The presence of system zeros or other 
poles near the root of interest makes the comparison to a 
second-order system invalid and reduces the accuracy of the 
damping ratio calculation. 
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Closely related to this method is a method based on a 
measured transfer function, calculated from cross- and power- 
spectra of the input and the measurement [14,15]. By analogy 
to the transfer function of a pure second-order system near 
the resonant peak, the damping ratio may be expressed as [9]: 





Im(T)co dco 


l o 



dco 


2 


where T = T(jto) is the measured transfer function. In 
practice, the transfer function often has other zeros and poles, 
therefore the limits of integration have to be limited about 
the natural frequency. If the damping ratio is small (0.01- 
0.05), the resonant peak frequency may be treated as the 
natural frequency. The integration limits 0.8 co^ to 1.2 co^ 
are often adequate, and can be further reduced if necessary 
to lessen the effects of nearby roots or noise [1]; with 
100 spectral lines in the range 0.8 co^ to 1.2 co^ , good 
damping values can be obtained by integrating from 0.95 co^ 
to 1.05 ojp. With 10-15 spectral lines between 0.8 co^ to 1.2(0^, 
the integration range 0.90 co^ to 1.10 co^ is adequate. Damp- 
ing ratio is thus given approximately by: 


? = 



9ce. 
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The principal effect of nearby modes is to shift the phase 
angle of the measured tranfer function. This effect can be 
minimized if, before computing £ by the above formula, the 
phase of the measured transfer function is changed by setting 



where T is the measured transfer function at the resonant 
P 

peak. 

Other methods of identifying transfer functions attempt 
to determine numerator and denomerator polynomials in jco 
by matching the experimentally determined magnitude and phase 
data [16]. The factored denominator then shows the system 
modes, and the damping of the mode of interest may be calcu- 
lated. These methods require specification of overall model 
order, however, and often fail to converge if the data contains 
noise. 

2.3 ADVANCED SYSTEM IDENTIFICATION TECHNIQUES 

System Identification is a broad technology which may be 
used to verify system simulations or develop new mathematical 
models of static or dynamic systems based on experimental 
data. These techniques are based on sound mathematical and 
statistical principles such that the resulting models and the 
estimated model parameter possess desirable statistical charac- 
teristics. The classical techniques based on heuristic ap- 
proaches are, therefore, being slowly replaced by more system- 
atic system identification methods. 

The overall system identification approach Is divided 
into pretest evaluation and algorithm verification and post-test 
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data processing. The post-test data processing consists of 
data conditioning, model structure estimation, parameter esti- 
mation and model verification (see Figure 2.1). The data 
conditioning involves filtering and removing undesirable 
characterist ics from data (e.g., spikes, outliers). The 
model order and the important linear and nonlinear terms are 
estimated in the model structure determination phase. The 
parameters of the selected model are identified in the next 
stage. The model verification stage applies several statisti- 
cal techniques to validate the overall model. A detailed 
description of the system identification process is given in 
Ref. 19. Other parameter identification methods are discussed 
in Refs. 20-27. 

Model order determination techniques select the simplest 
system which can explain the measured response adequately. 

The criteria for determining model adequacy include: (a) pre- 

diction error, (b) Akaike information criteria, (c) fit error, 
(d) F-ratio and several others. A detailed description is 
given in Refs. 30 and 31. 

The parameter estimation methods may be divided into two 
broad classes, the off-line approaches and on-line approaches. 
Several methods have been developed for off-line application 
but the maximum lieklihood method has been the most successful. 
The likelihood function of the output with respect to the 
parameters is maximized by selecting the parameter values. 

For numerical purposes, the negative logarithm of the likeli- 
hood function is minimized. The method optimally treats both 
the state disturbance noise and the measurement noise. Ac- 
curate estimates of all required parameters are obtained at 
the cost of high computation time and storage requirements. 

One of the more successful on-line parameter identifica- 
tion techniques is the instrumental variables approach. It 
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Figure 2.1 General System Identification Approach 



is an extension of classical least squares techniques, but 
accounts for all noise sources. It is discussed in more 
detail in the next chapter. 

System identification techniques determine all unknown 
parameters of the system directly and can therefore estimate 
frequencies and dampings of all modes. Therefore, modes that 
are close to each other in frequency may be isolated and the 
effects of noise in the process and the measurements are mini- 
mized. 

2 . 4 SUMMARY 

This chapter reviewed methods which have been used or 
attempted in the on-line determination of the stability 
of dynamic systems. Some are capable of being automated, 
while others require substantial subjective interpretation 
by the analyst. All except parameter identification 
methods place considerable weight on analogies to classical, 
well-understood, second-order dynamics, and are hampered by 
experimental conditions which obscure this analogy, such as 
noise or the presence of nearby modes. Parameter identifica- 
tion methods, on the other hand, can resolve these problems 
effectively, but they must be specialized to the problem 
of on-line dynamic stability measurement. 

The next chapter describes a nexv on-line parameter iden- 
tification algorithm, based on the instrumental variables 
method, that overcomes these problems. Its derivation and 
initial evaluation are discussed in the remainder of this 
report. 
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CHAPTER III 

A NEW REAL TIME TECHNIQUE FOR MODEL IDENTIFICATION 

3.1 INTRODUCTION 

From the previous discussion, it is clear that current 
approaches have several deficiencies. To obtain reliable esti- 
mates of frequencies and dampings of various modes from a 
limited wind tunnel test time, it is desirable to have an 
algorithm which can extract the maximum information from the 
data. Such an algorithm must be able to handle multi-input, 
multi-output responses, must correct for measurement errors 
(errors in instrumentation) and process noise errors (random 
wind tunnel turbulence) , and must consider interaction between 
many modes. An on-line or real time technique is required 
and it should be able to handle modes with low damping ratios 
(0.01 to 0.10), characteristic of structural modes near regions 
of aeroelastic instability. 

3.2 PROBLEM DISCUSSION AND MODEL REPRESENTATION 

The dynamic relationships between the inputs and the out- 
puts may be expressed in several forms. For the purpose of 
measuring dynamic stability in the 40- by 80-foot wind tunnel, 
a general model is most suitable since the mode shapes are not 
known a priori and may change with test condition. The model 
may be either: (a) modal state variable or autoregressive 

moving average (Appendix A) , (b) discrete time or continuous 
time, and (c) in time domain or in frequency domain. The ARMA 
representation is better than the state variable form because 
approximate values of natural frequencies and dampings 
(required in parameter estimation in model state variable rep- 
resentations) of the modes are not available. The continuous 
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time representation is used, in spite of discrete measurements, 
mainly for improved numerical conditioning, when the damping 
ratios are small (less than 0.10). The algorithm is developed 
in time domain and is subsequently converted into frequency 
domain for covenience of implementation. 

For the multi-input/multi-output systems, there are sever- 
al forms of the continuous ARMA representation (see Hannan [29] 
for corresponding discrete ARMA representations). We use the 
following 


2 a.I y Cn ' l} (t) 
i=0 1 


n l n 2 

2 B.u^ n ' x) (t)+ Z C.v^ n_1 ^(t) 

i=0 1 i=0 1 


0 < t < T 


a 


o 


1 , 



(3.1) 


where ysR p , ueR q , vsR p , B.sR pxq , C jL sR pxp , a-eR and (»)^ 
denotes the ith differential of (•)• The order of the model 
n is twice the number of modes, n^ is the number of zeros 
of the transfer function between the outputs and the inputs and 
n 2 is the number of zeros of the transfer function between the 
outputs and noise, n^ and n^ are usually equal to or less 
than n. Since a^'s are scalars, the poles of the system 
are roots of the equation 


n 

Z a. s 1 = 0 (3.2) 

i=0 1 

which depends only on parameters (that is why we use 

the form of Eq. (3.1)). Matrices together with a^ deter- 

mine the input/output transfer functions, i.e. 
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7 11 n - * 

y(s) = S B.s 11 ’ 1 u(s) 

7 1 A(s) i=Q i 

where 

A(s) = s 11 + a ,s n_1 + a 7 s n ~ 2 + ... + a„ 
l l n. 

3.3 PARAMETER ESTIMATION 

As mentioned in the previous section, the input/output 
transfer function depends only on coefficients and matrices 

B^. Therefore, in the measurement of dynamic stability it is 
of major interest to estimate these parameters accurately. 
Matrices give the noise-output transfer function and need 

not be estimated. The presence of correlated noise terms, 
however, precludes straightforward application of the least 
squares method. 

Instrumental variables is an ideal technique for such 
application. Developed initially by econometricians, it was 
studied in detail by Wong and Polak [23] and recently applied 
to vehicle estimation problems by Gupta and Hall [28]. To 
explain the technique, let us consider a static system, where 
the p outputs y are related to the pxm independent 
variables X and m parameters 0 by the equation 

y=XG+s (3.5) 

where s is a zero mean white noise with unit covariance. If 
e and X are uncorrelated , an efficient estimate (one with 
the smallest mean square error) is 

© LS = (X T X) _1 X T y (3.6) 


(3.3) 

(3.4) 
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are 


The mean and covariance of 0 

E(e LS ) = e 

E[E(0 ls ) - 9] [E(0 ls ) - e] = (X T X) (3.7) 

where E(*) is the expected value of (•)• However, if e and 
X are correlated, the least square estimator is biased 

B(0 ls ) = 0 + (X T X)' 1 E (X T s ) (3.8) 

An instrumental variables estimator could give unbiased esti- 
mates. In this method, we select a pxm matrix, Z, which 
is uncorrelated with s. Then, 

§ IV = (Z T X) Z T y (3.9) 


This estimator is unbiased. The accuracy of the estimate 
depends upon the choice of Z. 

The application of the instrumental variables technique 
to the estimation of ARMA parameters is described in Appendix 
B and shown schematically in Figure 3.1. The following high- 
lights of the algorithm are especially pertinent: 

(1) The estimation of parameters using the instrumental 
variables approach requires the solution of a set 
of linear equations. These equations have a very 
specific form for the multi-input, multi-output 
case. This specific structure must be utilized to 
save computation time and storage. 

(2) The elements of the matrix in the linear equations 
depend upon the Fourier transforms of inputs and 
outputs. A fast Fourier transform should be used 

for this purpose and, if feasible, 2° data points 
should be collected, where a is an integer. 
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INPUT & OUTPUT 
TIME HISTORIES 


F REQUENCY RANGE OF INTEREST 
TIME DELAY PARAMETER 


NO. OF ZEROS AND NO. OF POLES 


FAST FOURIER 
TRANSFORM 


SET UP 

SIMULTANEOUS 
LINEAR EQUATIONS 


FOURIER TRANSFORM OF 
INPUTS & OUTPUTS 


A, b 


▼ 

SOLVE LINEAR 
EQUATIONS 


TRANSFER FUNCTION 
NUMERATOR AND DENOMINATOR 
POLYNOMIALS 



Figure 3.1 Schematic Flow Chart of the On-Line Dynamic 
Stability Measurement Algorithm 
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(3) The algorithm is flexible in that the frequency range 
of interest may be selected a priori. 

(4) The inputs and outputs may be scaled for better 
numerical conditioning without affecting pole and 
zero locations. 

3.4 MODEL ORDER DETERMINATION 

In real time estimation of dynamic stability, the number 
of important modes is not known, particularly in complex aero- 
elastic models. Therefore, it is often necessary to estimate 
the number of significant modes in the output measurements. 
Methods to determine model order have been described by 
Akaike [30] and Tse and Weinert [31], These techniques are 
applicable for systems driven by unknown random inputs. 

We propose here a fundamentally different approach to the 
problem of estimating the order of the model based on input/ 
output measurements. The basic idea is as follows. If the 
specified model order is lower than the true model order, the 
identification algorithm in this chapter will give incorrect 
estimates of poles and zeros. If the specified model order 
is the same as the true model order, correct poles and zeros 
will be identified. However, if the specified model order is 
higher than the true model order, there will be many superflu- 
ous poles and zeros in addition to the true poles and zeros. 

The superfluous poles must approximately cancel out the super- 
fluous zeros. Therefore, if we specify a maximal model order 
in the algorithm of this chapter, the true model order could 
be obtained after the transfer function is simplified by 
cancelling the poles and zeros which are very close to each 
other. 

The identified transfer function may be simplified as fol- 
lows. First, poles and zeros which are very close to each other 
may be eliminated reducing the order of the transfer function. 


20 



Second, the poles and zeros outside the frequency range of in- 
terest may be removed. Let the identified transfer function be 



3.5 APPLICATION OF THE MODAL ESTIMATION ALGORITHM TO GROUND 
RESONANCE SIMULATION DATA 

A simplified model for the ground -resonance phenomenon 
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(3.12) 


where 

£ lc , ? ls are cyclic lag modes; 
y is lateral hub motion; 

F is the shaker input; 

Ij., , lag inertias = 1; 

My, support mass = 30 
Cy, support damping = 0.05 
ay, support frequency = 5 Hz 
, lag frequency = 0.25 
, lag damping; and 
fl, the rotor speed = 10 Hz. 


Time has been nondimensionalized with the rotor speed.. A linear 
simulation was conducted at 10 Hz rotor speed with the lag damp- 
ing of 31.4. Using Eqs . (3.12), the poles and zeros of the trans- 
fer function between the lag modes and lateral motion vs. the 
shaker input are shown in Table 3.1. A sine sweep input is 
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Table 3.1 


Poles and Zeros of the Ground Resonance Simulation 


POLES 


MODE 

EIGENVALUE 

NATURAL 

FREQUENCY 

DAMPING 

FACTOR 

Upper-Lag 

-.237 + 1.164j 

1.20 

.239 

Lower-Lag 

-.224+ • 86 7 j 

.896 

.250 

Lateral Hub 

-.0230 + . 503 j 

.503 

.0457 



OUTPUT C ]c 

OUTPUT C- 

LATERAL HUB POSITION, Y 


-1.25 

-.25 

-.271 +7.00J 

ZEROES 

.75 

0.0 

-.229 + l.OOg 


0.0 

0.0 



0.0 



GAIN 

.0345 

.069 

.0345 


applied and 512 measurements are taken at intervals of 0.5. 

The algorithm of Section 3.3 is applied to determine the effects 
of: (a} time delay parameter t in the instrumental variables 

approach, (b) output noise level, (c) outputs used in estima- 
tion, (d) number of data points, and (e) discrete frequency 
noise, on poles and zeros estimation accuracy. In the applica- 
tion of the algorithm, the nondimens ional frequency range of 
interest is limited to 0.4 < co < 1.2 and the numbers of poles 
and zeros in each transfer function are specified to be six 
and five, respectively. 
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Effect of Time Delay Parameter t 


Table 3.2 shows the estimated poles and zeros between 
outputs ^ an ^ y and “the shaker input as a function 

of the time delay parameter t when each output is contamin- 
ated by 101 random noise. Since the system is of order 6, the 
instrumental variables approach requires a value of t exceed- 
ing 12 sample points or 6 seconds. The estimates of the poles 
are shown in Figure 3.2. Note that the estimates of the poles 
are most accurate at t = 10 sec. As explained before, the 
estimated transfer functions for r = 10 may be simplified to: 


^lc^ S ) . 0432s 2 (s - . 515) 
F ( s ) A 


gls (s) , . 0117s(s + . 589) (s .667+.205j) 
F (s ) A 


(3.13) 


y (s) _ . 0 29 1 (s + .245 t . 937j)(s + . 202 ± 1.05j) 
F (s) A 


A = (s + .0233 + .50 3j) (s + .221+.854j)(s+ . 2 67 + 1. 15 j) 

(3.14) 

Note that the estimates of poles are more accurate than those 
of the zeros. This is because all three outputs contain in- 
formation regarding the poles while the information about the 
zeros of any particular transfer function is contained only in 
one output. 


Effect of Output Noise Level 

Keeping the time delay parameter fixed at 10, the output 
noise level is varied. Table 3.3 shows the poles and zeros 
for various noise levels up to 100% (see also Figure 3.3). 
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Table 3.2 


Effect of Time Delay Parameter t on Estimates of 

Poles and Zeros 



SIMULATION 

T = 0 

x = 2 

T = 10 

x - 20 



-.287 + I.164J 

-.124 + 1.098J 

-.149 + 1 . 13j 

- .267 + 1.1 5 j 

-.285 + T.124J 

POLES | 

-.224 i . E67j 

-.122 + . 875j 

-.751 + .871 j 

-.221 + ,864j 

-.223 + .854j 



-.0230 + - 503 j 

-.0245 ± . 504j 

-.0244 + ,504j 

-.0233 t -503j 

-.0232 ± .503j 

S is) 


-1.25 

-.218 + 1.02j 

-.334 + l.llj 

-1.64 + 1 .44j 

-.896 

Tc' ' 
F(s) 

ZEROES 

.75 

.285 + , 439 J 

.302 + .401 j 

.513 

+.698, 175.7 


0.0 9 Q.O 

-.591 

-.481 

- .0130, J14 

-.0157+ .0687J 


GAIN 

— 

.0359 

.0309 

.00907 

-.000191 




-.202 + .988j 

-.334 + 1.02J 

.667 + . 205 j 

.266 

ZEROES 

0.0 

.109 + . 333 j 

.108 + .361 j 1 

-.589, -.0387 

-.191 + .280j 

■mi 



4.0 ' j 

4.56 

7.47 

-2.93, 2.83 

■ 

GAIN 

.069 

-.0113 j 

-.00847 

-.00156 

-.00605 

ESI 

1 

ZEROES 


-.0867 + 1 .Q9j 
-.103 + . 882 j 

-.112 + l.lOj 
-.132 + .889j 

-.202 + 1.05j 
-.245 t 937 j 

-.18 + K04j 
-.287 + . 891 j 




21.0 

20.2 

15.0 

11.8 

i 

GAIN 


-.00147 

-.00146 

-.00194 

-.0023 



Figure 3.2 Effect of t on Estimates of Poles (Complex 
Conjugate Poles Not Shown) 
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Table 3.3 

Estimates of Poles and Zeros for Various Output Noise Levels 





























































Many of the estimated transfer functions may be simplified as 
shown before. 

Effect of Outputs Used in Estimation 

For a 10% noise in outputs and t =10, Table 3.4 shows 
the effect of outputs on the estimation accuracy of poles and 
zeros. Note that none of the outputs can individually estimate 
all the poles accurately. However, a combination of outputs, 
particularly a-nd Y> gives excellent estimates of 

all poles and a reasonable estimate of the zeros. 

Effect of Number of Data Points 

The estimates of poles and zeros obtained using 512 points 
are compared with the estimates obtained using only 128 and 
256 points in Table 3.5. The sampling interval is kept fixed 
at 0.5, the noise level is 10% and t is 10. The estimate of 
the damping factor for the lateral hub mode deteriorates quite 
significantly when the estimate is based only on 128 points. 

Effect of Discrete Frequency Noise 

Table 3.6 gives the estimates of poles and zeros when the 
outputs have various levels of one per rev noise in addition 
to the 10% random noise. Again r is set at 10 and the fre- 
quency range of interest is set at 0.4 < co < 1.2. A 10% one 
per rev noise causes a large error in the estimates of the fre- 
quencies and dampings of the system modes. 

Next, the frequency range of interest is reduced to 
0.4 < co< 0.95 and 1.05 < co< 1.20 to rej ect the one per 
rev noise. The new estimates are shown in Table 3.7. These 
estimates are more accurate than the estimates of Table 3.6. 



Table 3.4 

Estimates of Poles and Zeroes of Transfer Functions for Varying Sets of 

Outputs Used in Estimation 
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Table 3.5 


Effect of Length of Data Estimates of Poles and Zeros of 

Transfer Function 



SIMULATION 

123 POINTS 

256 POINTS 

512 POINTS 



-.0230 +- 503 j 

-.0389+ .496 j 

-.0232+ .503j 

-.0233+ .503 j 

POLES 

-.224 + .8S7j 

-.223+ .855 j 

-.235+ .862 j 

-.221 + .864j 



- . 287 + 1.1 64 j 

-. 308 + 1. T3j 

-.307+lJ3j 

-.267 + 1.15J 



-1.25 

- .0445 

-.185 

-.0130 

ZEROES 

.75 

-.382+ .337j 

.238 

.114 

PCs) 


0.0 

.967 

.507 

.513 



0.0 

4.95 

-23.3 

-1.64 + 1. 44j 


GAIN 

.0345 

--■0119 

.00241 

.00907 

q s (s) 

ZEROES 

-.25 

0.0 

.314 

-.208+ .433 j 

.264 

-.200+ .273 

.667+ .205 j 
-.589 

FCs) 


0.0 

3.97' 

-2.91 

-.0387 




4.27 

3.15 

7.47 


GAIN 

.069 

-.00562 

-.00877 

-.00156 

xisl 

F(s) 

ZEROES 

-.Z71+_1.00j 
229 + 1. 00j 

-.Z18 + .8?9j 
-.118 + 1.08j 
1.78 

-.297+ .830J 
-.135 + 1.05j 
4.32 

-.245+ .937j 
-.202+ 1.05J 
75.0 

' 

GAIN 

.0345 

-.0303 

-.0124 

-.00194 


Table 3.6 

Effect of 1/Rev Discrete Frequency Noise on Estimates of 
Poles and Zeros (Random Noise 10%, r=10, . 4 < to < 1.2) 


• 

SIMULATION 

502 1/REV 
NOISE 

20% 1/REV 
NOISE 

IG£ 1/REV 
NOISE 

U 1/REV 
NOISE 

POLES 

-.0230 + .503j 
-.224 + .867J 
-.287 + 1.164j 

-.0350 + .497j 
-.0621 + T.OOj 
-.868 
1.17 * 

.113 

-.0286 + . 497 j 
-.144 + .988j 
-.693 

-.0232 + .502j 
-.416 + .726j 
-.213 + .993j 

-.0233 + . 503j* 
-.227 + .864j 
-.278 + 1.T4J 

C lc (s) 

ZEROES 

-1.25 

.75 

0.0 

0.0 

.171 + .591 j 
-.307 + . 954j 
T .79 

.0748 

-.846 + . 223 j 
.259 + .593j 

-.314 + .324j 
.310 + ,422j 
-7.91 

-.0120 

.110 

.520 

-2.16 + 1.06J 

n 

GAIN 

.0345 

.0817 

.0490 

.0082$ 

.00727 

? ls (sl 

ZEROES 

-.25 

0.0 

G.O 

-.300 

.204 + .619 j 
-.348 + .987 j 

.211 + .492j 
-.318 + . 4 T 7j’ 
2.94 

-.!76+.Z76j 

.336 

-.759 

1.31 

.531 + .1743 
-.473 
-.354 
5.12 


GAIN 

.069 

.148 

-.0281 

-.0315 

-.00264 


ZEROES 

-.271 + T.OOj 
-.229 + T.OOj 

.0514 + .918j 
-.849 

.986 + 1.58J 

.217 

-.0514 + .942j 

-.942 

-77.16 

-.550 + .689 j 
-.101 + .956j 
-37.8 

-.269 + .963j 
-.198 + l.OTj 
16.17 

GAIN 

.0345 

-.0165 

.000385 

.000733 

-.00180 
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Table 3.7 


Effect of 1/Rev Discrete Frequency Noise on Estimates of 
Poles and Zeros (Random Noise 10%, t =10, 

.4 < to < . 95, 1.05 < to < 1.20) 




SIMULATION 

50% I /REV 

20% 1/REV 

10% 1/REV 

1% 1/REV 



NOISE 

NOISE 

NOISE 

NOISE 



-.0230 + . 503j 

-.0223 + .504j 

-.0232 + . 503j 

-.0234 + . 503j 

-.0235 + .503j 

POLES 

-.224 + .867j 

-.182 + .660j 

-.240 + .802j 

-.223 + . 834 j 

-.207 + . 850j 



-.287 + 1.164J 

-.184 + 1.02j 

-.238 + 1.06J 

-.247 + l.lOj 

-.237 + 1.1 4 j 



-1.25 

-.130 + ,542j 

.356 + .281 j 

.354 + . 1 79 j 

-.191 

'let*) 

ZEROES 

.75 

.297 + . 506j 

-.359 + . 388j 

-.351 

.344 + . 132j 

■iQK 


0.0 

-3.28 

-6.75 

-.863 

-1.29 + 1 .21 j 



0.0 



-3.56 


■I 

GAIN 

.0345 

.0197 

.00654 

.00831 

.0120 



-.25 

.150 + .344j 

-.0516 

.0463 

.0883 + .179j 

ZEROES 

0.0 

-.167 + .618j 

.275 

.131 

-.953 + .495j 

tt?t 


0.0 

1.65 

-.444 + . 51 9 j 

-.666 + . 474 j 

3.28 




1.46 

2.00 




GAIN 

.069 

-.0373 

-.0253 

-.0153 

-.00706 

mm 


-.271 + 1.00J 

-.185 + * 605 j 

-.358 + 784 j 

-.347 + .91 0J 

-.172 t *940j 

11 

ZEROES 

-.229+l.OOj 

-.0183 + . 940j 

-.0768 + . 967 j 

-.115 + .973j 

-.254 + 1.06J 



5.43 

7.22 

14.98 

-246.8 


GAIN 

.0345 

-.0800 ! 

-.00452 

-.00193 

.000109 


3.6 APPLICATION OF THE MODAL ESTIMATION ALGORITHM TO SIMULA- 
TION DATA OF A CANTILEVER WING WITH ROTORS AT EACH WING 
TIP 

The symmetric motions of a cantilever wing with rotors 
at each wingtip are approximately described by a 9-DOF model 
discussed by Johnson [5]. Of the four possible inputs, xve will 

restrict our discussion to the collective pitch input, 9 . 

Two outputs which will be considered are the vertical wing 
bending, q^, and chordwise wing bending, q 2 - The poles of 
the system and the zeros between the outputs q^ and q£ and 
the collective pitch input are shown in Table 3.8. Note that 
in the q./0 Q transfer function, four complex zero pairs 
approximately cancel four pole pairs, while in the q 7 /0 

L O 
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Table 3.8 

Symmetric Modes of a Cantilever Wing 


MODE 

POLE 

NATURAL 

FREQUENCY* 

DAMP TNG 
FACTOR , C 

Coning (S 0 ) 

-0.2095 +2.55<l3j 

2.5533 

.02174 

Upper Inplane (4^) 

-.07232 + 2.4345j 

2.4355 

.0299 

Upper Out-of-Plane 

(e H ) 

-.3015 i 1.8546 j 

7.3789 

.1605 

Wing Torsion (p) 

-.07785 i 1.344J 

1.3452 

.0573 

Wing Cncrdwise 

V 

Wine Vertical 

-.02339 + .5557j . 
-.01914 +_ . 3935 j 

.6667 

.399 

.0426 

.04798 

Lower Enplane 

fc-i> 

1 

- -.03348 + .-23563 

.3471 

.255 

Lower Out-of-Plane 

CB.,) 

-.372 + * 1 763j 

.3534 

.871 

? o 

-.7445 



h 

-.04379 





OUTPUT q 1 

OUTPUT q a 


-.0753 t 2-435J 

-.32? £ 8.73j 


.0111 £ 2. 05j 

-.0673 £ 2. 43 j 


-.297 £ 1.338 j 

-.303 £ 1.85j 

ZEROES 

-.147 + 1.063j 

-.0750 £l.34j 


.0762 £ .627 j 

-.311 £ - 1 73 J 


-.340 £ . 145j 

-.0195 £ .3 99j 


-.0939 £ . 323 j 

-.0890 £ - 338j 


-.0759 

-.00305 


0.0 

0.0 

GAIK 

-.0727 

.00198 


*Per rev 

transfer function, six zero pairs approximately cancel six 
pole pairs. In fact, in the frequency range of interest, 
0.2 < w < 0.75, the transfer functions c l 1 / e 0 anci ^ Q o 
may be written as (see Eqs . (3.10) and (3.11)) 

q l t ' S ' ) _ _ . 00498 (s - . 0762 + . 621) ) 

0^ (s) (s + . 0284 + .666j)(s + .0191 • 3 9 8 5 j ) 
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q 2 (s) _ -0230 

0 q (s) (s + .0284s + .666j) 

A simulation was conducted using a sine sweep collective 
pitch input. The outputs q-^ and q 2 and the input 0 Q were 
sampled at intervals of 0.5. 512 samples are used in the iden- 

tification. The frequency range of interest is 0.2 < w < 0.75 
and a 5-DOF model is specified. The identified poles and zeros 
for two cases are given in Table 3.9. It is clear that many 
estimated poles are close to estimated zeros. When output 


Table 3.9 

Identified Poles and Zeros for Cantilever 
Wing Simulation Data 

5 Modes, t =10, 512 Points, 0.5 Sec Sampling 
0. 2 < w < 0. 75 





OUTPUTS USED 

IN ESTIMATION 



SIMULATION 

q 2 

AN0 



-.0234 + . 666 j 

-.0257 + . 566 j 

-.0233 * .661 j 



-.0191 + . 398j 

.0066 + . 437j 

-.0137 + . 395 j 

POLES 

-.0885 + . 335j 





-.312 + .I76j 
-.144 

-.036 + . 1 55 j 

• 

-.011 + .129j 



-.0432 

0.0 + .694j 

-.002 + .691 j 




•Oil + . 6 1 6 j 

•015 i .5320 



-.0196 * . 399j 

-2.72 

-2.31 



-.0311 * J73J 

0.0 + .691 j 

-.005 + . 537 j 

q 2 (s) 

ZEROES 

-.0889 + . 336 j 
0.0 

.010 + .615j 

.0119 + - 532 j 



0.0 

.0028 + . 435j 

-.0193 1 . 393 j 




0.0 + . 1 55 j 

.0318 + . 1 54j 


GAIN 


.00793 

.00799 


| 

! .0762 + .S21j 


I -2.06 

q-j(s) 


| -.340 + . 1 45 j 


-.0037 + ,592j 

| ZEROES 

j 

-.0939 + . 323 j 
-.0159 


.0671 + , 623 j 

.0154 t • 53 6 j 


1 

0.0 


.0382 + . 1 50 j 


GAIN 



-.00206 
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c^2 is used In estimation, the estimated transfer function is 
approximately 

q 2^ = .0216 

0 (s) (s ■*■ . 025 7 +_ .666j) 

When the outputs q^ and q 7 are used in estimation, the 
transfer functions may be simplified to 

q l ! - S ^ _ .00424(s - .0671 + .623j) 

6 (s) “ (s + . 0283 + . 661j ) (s + . 0137 + . 395j) 

q 2^ S ^ _ .0225 

0 (s) (s + .0283 + . 661 j ) 

This shows that when the specified model order is much higher 
than necessary, enough poles and zeros cancel out to give a 
low order transfer function. Also, the estimate of the chord- 
wise wing bending mode improves when both outputs q^ and 
are used in estimation. 

3.7 CONCLUSIONS 

A new algorithm based on the instrumental variables ap- 
proach is developed for dynamic stability estimation from 
input/output signals. A preliminary evaluation of the algo- 
rithm on simulation data of low order models demonstrates its 
effectiveness in estimating frequencies and damping factors of 
important modes. 
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CHAPTER IV 

APPLICATION TO XV- 15 SIMULATION MODEL 


4.1 INRODUCTION 

The on-line modal identification algorithm of Chapter III 
was demonstrated on low-order models. This chapter describes 
the application of this algorithm to high- order simulated data 
from an XV- 15 math model and examines both algorithm character- 
istics per se , as applied to this model, and its ability to 
identify XV-15 aeroelastic modes of interest across the air- 
craft's flight envelope. 

The following sections present a description of the XV-15 
math model used in this investigation, and discuss the ability 
of the algorithm to identify prominent XV-15 structural modes 
from measurements of the model states or, alternatively, from 
measurements (obtained from a typical wind tunnel test) com- 
prised of one or several model states. The math model is con- 
sidered in Section 4.2. The application of the algorithm to 
time histories of symmetric and antisymmetric motions is 
described in Sections 4.3 and 4.4., respectively, and a sum- 
mary of results is given in Section 4.5. 

4.2 XV-15 ANALYTIC MODEL 

The XV-15 model used in this study is a linear simulation 
of the elastic airplane mounted on the model support system of 
the NASA full-scale (40- x 80-ft) wind tunnel. The model, 
furnished by NASA and described in Ref. 5, utilizes predicted 
XV-15 structural characteristics and configuration parameters, 
and contains a total of 13 degrees-of-freedcm: ten of the 

aircraft and three of the support system, in each of a set of 
symmetric equations and a set of antisymmetric equations. 


fBECSMNG PAGE BLANK NOT FILMED 
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These degrees of freedom, and their associated modes, are 
identified in Appendix C, Table C.l. Since each equation set 
contains 13 modes, it has a characteristic equation of 26th 
order. This model has produced, in previous investigations, 
results showing very good agreement with experimental data 
[3-5]. 

Typical XV- 15 dynamic characteristics are shown by the 
symmetric roots plotted in Figure 4.1. A number of low-damped 
oscillatory roots are seen along the jco-axis; these roots 
include the wing bending and torsion and wind tunnel balance 
modes. Higher-damped rotor modes are seen farther left of the 
jco-axis. Typical damping ratios of the former are on the order 
of 1 to 2 %. These roots change with airspeed and airplane 
configuration (pylon angle) , and as illustrated by arrows in 
this figure, tend to move toward the jco-axis with increasing 
speed. 

NASA-derived dynamic model data used here consists of 
the coefficient matrices of the dynamical equation, in NASA 
formulation, 

Ay'x + A-^x + A q x = Bu (4.1) 

where x and u are described in Table C.2 (Appendix C) , the 
A o> A-^, and A 2 matrices contain system aerodynamic, struc- 
tural, and inertial characteristics , and B is a control or 
disturbance input effectiveness matrix. The elements of these 
matrices vary according to flight condition. 

This equation may be easily rewritten to treat the degrees 
of freedom as state variables: 
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Figure 4.1 Typical Characteristic Root Locations, 13 Degree- 
of -Freedom Model; XV- 15 Plus Wind 
Tunnel Supports 


x = Fx + Gu + Tg 



where x, y, u, and g are defined in Table C.2 (Appendix 
C) , and F, G, F, and H are the system, control, gust 
disturbance, and observation matrices, respectively. From this 
analytic formulation are computed the system eigenvalues, 



transfer functions, and time histories used in this study. 
The analytic results are, of course, identical to those from 
the NASA equation formulation. 


4.3 APPLICATION TO XV- 15 SIMULATION MODEL: SYMMETRIC MODES 

Six flight conditions representative of XV- 15 free-air or 
wind tunnel testing are selected as test points at which the 
XV-15's dynamic characteristics would be computed and the 
algorithm's effectiveness in identifying roots investigated. 

The six flight conditions are shown in Table 4.1. This exercise 
applies the algorithm to time history data containing thirteen 
dynamic modes and is representative of the application for 
which the algorithm was developed. 

The approach taken here in generating and identifying 
XV-1S dynamic characteristics is illustrated in Figure 4.2, 
indicating that simulation time histories are computed as a 
substitute for experimental data. The objective of this 
exercise is twofold: (1) to investigate the basic character- 

istics of the algorithm per se and the effects of algorithm 


Table 4.1 

Selected XV-15 Flight Conditions 


CONDITION 

V KTS 

Op(°) 

a, rpm 

1 

0 

90 

565 

2 

80 

82.3 

. 565 

3 

100 

60 

565 

4 

120 

30 

565 

5 

160 

0 

458 

6 

200 

0 

458 
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Figure 4.2 XV- 15 Simulation Mode Identification Procedure 


modifications on mode identifications, when applied to high- 
order data; and C2) to investigate the algorithm's ability 
to identify specific XV-15 modes that would be of interest in 
a test situation. In each case, the results of interest are 
the accuracy of identification of frequency and damping of 
critical modes and the overall placement of transfer function 
poles and zeros, judged by comparison to the analytical roots 
determined from the state equations. In the first part, 
measurements of airplane states are used directly for identi- 
fication; in the second, the time histories of horizontal 
and vertical wingtip deflections are computed and used as 
measurements, more realistically simulating experimentally 
available data. 

To evaluate basic algorithm characteristics, the simula- 
tion model at flight condition 3 is excited by a random 
collective pitch input and the resulting time histories of 
wing vertical and chordwise bending motions (q^ and q£, 
respectively) are used as measurements. The random control 
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input is white noise, filtered with a low pass Martin-Graham 
filter having a constant response up to 0.94 rad/sec and 
1.57 rad/sec cutoff frequency. The measured input frequency 
spectrum is shown in Figure 4.3; the corresponding autocor- 
relation characteristics of the input time history are shown 
in Figure 4.4; and an example of this time history is shown 
in Figure 4.5. 


It should be noted that all frequencies in this study are 
non- dimensionali zed with respect to rotor RPM; thus, a fre- 
quency of 1 rad/sec actually represents a frequency of 

rad/sec , 


Q, (RPM) u _ __ 2nQ 

s' A n 


60 


or 


60 


true scale. 


The effects of the following algorithm parameters are 
investigated: (1) assumed model order; (2) time delay para- 

meter, t ; (3) specified number of zeros; (4) the use of two 
measurements; (5) the effect of specified frequency range; and 
(6) the effect of various forms of scaling within the algorithm. 
The results obtained are as follows: 


(1) Model Order. Altering the number of poles and zeros 
specified for the algorithm (i.e., altering the 
number of degrees -of - freedom) changes the locations 
of all the roots. Figure 4.6 shows the effect of 
assumed model order on the damping of modes identified 
from q^ and q ^ measurement time histories. Best 

identification of this system was obtained with modal 
orders of 8 to 12. 

(2) Time delay parameter. The recommended value of t 
• depends on model order and sampling interval, as 

t = 2*n*At . 

Figure 4.7 shows that using values smaller than 
recommended may lead to large errors in damping. 

(3) Number of zeros. The number of zeros specified is 
generally equal to, or one less than, the number of 
poles. The damping values are sensitive to this 
parameter, as seen in Figure 4.8, and different modes 
may be affected in different ways . 
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AMPLITUDE 


FILTER 


NOTES: • FIRST 40 SEC OF 256-SEC RECORD 

• FILTERED WHITE NOISE 



Figure 4.5 Example of Random Control Input Time History 


(4) Frequency range. Restricting the frequency range 
over which the algorithm identifies a model may im- 
prove the identification by eliminating extraneous 
frequencies, but it also reduces the amount of infor- 
mation available. The specified range is thus a 
compromise; in general, a wide range (e.g., 

0.1 < to < 0.9) is desirable, and frequency scaling 
in some instances yields improved results. This is 
discussed next. 

(5) Scaling. Because of several numerical operations 
performed in the algorithm, the magnitudes of the 
quantities considered can affect identification 
accuracy due to computer roundoff and program trunca- 
tions. Specific conclusions are: 

(a) Scaling FFT output within the algorithm by a 

large multiplicative constant may be required to 
ensure convergence of external equation-solving 
and rooting routines used by the algorithms 
implying the need to "tailor" the algorithm to 
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Figure 4.8 
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• t = 10 
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Effect of Specified Number of Zeros on Predicted 
Damping; 13 DOF Model 
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the computer system it is run on. Further, 
scaling the measurement and control input FFT's 
by their respective RMS values was found to 
lead to a substantial increase in accuracy in 
the two -measurement case, as seen in Table 4.2. 

(b) Scaling the frequency range such that the maxi- 
mul value is unity leads to better mode identi- 
fication in some cases through computational 
improvement . 

(c) Scaling FFT output to give a higher weighting 
to data near the natural frqeuency of modes 
improves identification accuracy. This is il- 
lustrated by comparing the root predictions of 
Figure 4.9, without scaling, to those of Figure 
4.23 (which will be discussed below) with such 
scaling . 

Figures 4.10 through 4.17 show a series of trial identifi- 
cation runs using the q 2 measurement and varying the para- 
meters discussed in items (1), (2), and (4) above, illustrate 
that, in this case (where the poles are well separated from 
other poles or zeros) : (a) decreasing the frequency range 

improves the pole identifications, and (b) successive reduc- 
tions in assumed model order yield nearly equivalent pole 
position estimates until the sixth-order model (three degrees- 
of-freedom) is reached, at which point the roots become poorly 
identified. 


Table 4.2 

Effect of Measurement Scaling 


MEASUREMENT 

ORDER 

h 

q 2 

SL 

q 2 UNSCALED 

10 

-.0033+. 295 j 

-.0080+. 55 lj 

-.0083+. 392 j 

q 2 SCALED 


-.0035+. 295 j 

-.0092+. 59 lj 

-.0092+. 392 j 

q 2 + q 2 UNSCALED 

12 

-.0085+. 291 j 

— 

— 

qj + q 2 SCALED 


-.0095+. 294 j 

-.0031+, 549j 

-.0069+. 396 j 

ANALYTIC VALUE 


-.0097+. 292j 

-.0106+. 573j 

-.0098+. 397j 
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Figure 4.9 XV-15 Simulation Symmetric Mode Identification 

Without Magnitude Scaling 
(Compare to Fig. 4.22) 
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Figure 4.10 XV-15 Simulation Symmetric Mode Identification 

Based on q ? Measurement 
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Figure 4.11 XV- 15 Simulation Symmetric Mode Identification 
Based on q2 Measurement: Effect 

of Reduced Frequency Range 
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4.12 XV- 15 Simulation Symmetric Mode Identification 
Based on q 7 Measurement: Varying Model Order 
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Figure 4.13 XV-15 Simulation Symmetric Mode Identification 

Based on ^ Measurement: Varying Model Order 
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4.14 XV-15 Simulation Symmetric Mode Identification 

Based on Measurement: Varying Model Order 
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Figure 4,15 XV-15 Simulation Symmetric Mode Identification 

Based on Measurement: Varying Model Order 
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4.16 XV-15 Simulation Symmetric Mode Identification 

Based on Measurement: Varying Model Order 
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Figure 4.17 XV-15 Simulation Symmetric Mode Identification 

Based on q 2 Measurement: Varying Model Order 
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From these and a large number of additional trial identi- 
fication runs, it is concluded that the assumption of a large 
model order, e.g., 10 or 12, represents a good starting point 
for identifying modes of this system. Reductions in model 
order and selection of frequency range are somewhat dependent on 
the particular root locations of a given flight condition, and 
different simplifications may be obtained in different flight 
conditions. The application of the algorithm for greatest 
accuracy is thus empirical to the extent of requiring engineer- 
ing judgment in the selection of these algorithm parameters. 
Additional development effort is being directed toward removing 
the need for this operator interaction. 

The algorithm was next applied to simulated time histor- 
ies of symmetric wingtip deflection, of which horizontal and 
vertical components were computed from the values of state 
variables using the equations and constants described in 
Table C.3, Appendix C. The actual states, as defined in this 
study, cannot be measured directly in an experimental situa- 
tion, whereas wingtip deflections may be calculated using 
accelerometers or strain gage instrumentation on the wing. 

The model was excited by random inputs to either collective 
pitch or flaperon deflection, using the same random noise and 
filter as above. This is, again, representative of an experi- 
mental condition. Mode identifications were performed in each 
of the six flight conditions identified in Table 4.1. 

The symmetric XV- 15 modes of most interest during wind 
tunnel testing are the wing vertical and chordwise bending, 
wing torsion and longitudinal support strut modes. The root 
locations of these and other modes are plotted in Figures 4.18 
to 4.28, with system zeros shown for each of the collective 
pitch and flaperon input cases. The movements of poles and 
zeros due to changing flight conditions and configurations 
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Figure 4.18 XV-15 Simulation Symmetric Mode Identification 

Based on ^IP Measurement, 6 0 Input. Hover. 
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Figure 4.19 XV-15 Simulation Symmetric Mode Identification 

Based on X^p Measurement, 0 Q Input. Transition. 
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Figure 4.21 XV- 15 Simulation Symmetric Mode Identification 
Based on X^jp Measurement, 0 Q Input. Transition. 
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Figure 4.25 XV-15 
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Simulation Symmetric Mode Identification 
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Figure 4.26 XV-15 Simulation Symmetric Mode Identification 

Based on Z^p Measurement, Sp Input. Cruise. 
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Figure 4.27 XV-15 Simulation Symmetric Mode Identification 

Based on X,pjp Measurement, 0 Q Input. Cruise. 
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may be noted. Complete tabulations of the transfer function 
poles and zeros are shown in Appendix C, Tables C.4 to C.9. 

The roots identified by the algorithm from the Xyjp 
and Zjjp time histories are shown superimposed on the analy- 
tical roots in Figures 4.18 to 4.28. Again, only roots in the 
low frequency range of interest are shown. The algorithm 
contains a wide frequency range specification (0.1 < to < 0.9) 
and utilizes scaling as discussed in items 5(a) and 5(c) 
above. In nearly all cases, good accuracy is seen in the pole 
locations for the modes of greatest interest, while somewhat 
more scatter is seen in the zeros. 

4.4 XV- IS MODEL IDENTIFICATION: ANTISYMMETRIC MODES 

The roots of the antisymmetric modes at each of the six 
selected flight conditions are shown plotted in Figures 4.29 
to 4.39, for the frequency range of interest. Complete 
tabulations of transfer function poles and zeros are given in 
Appendix C, Tables C.10 to C.15. 

The antisymmetric roots of most interest are the wing 
vertical and chordwise bending, wing torsion, and strut lateral 
(side) deflection. As in the symmetric case, these roots tend 
to become less stable as airspeed increases. 

Figures 4.29 to 4.39 show the identification of these 
roots by the algorithm, superimposed as before on the analyti- 
cal roots. 

4 . 5 SUMMARY 

This chapter considers the ability of the algorithm dis- 
cussed in Chapter III to identify specific modes from time 
history data containing many degrees of freedom (in this case, 
13) . An analytical model of the XV-15 aircraft mounted in the 
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Figure 4.29 XV-15 Simulation Antisymmetric Mode Identifica- 

tion Based on Z T jp Measurement, Input. Hover. 
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Figure 4.30 XV-15 Simulation Antisymmetric Mode Identification 

Based on Xpjp Measurement, 0 Q Input. Transition. 
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Figure 4.31 XV-15 Simulation Antisymmetric Mode Identification 

Based on Z TIp Measurement, 8 p Input. Transition. 
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40- by 80-foot wind tunnel is used to provide representative 
simulated time history data. The conclusions from the work 
of this chapter are: 

(1) In most of the flight conditions examined, important 
aeroelastic modes are identified with good accuracy 
from time history measurements of wingtip deflection. 

(2) The algorithm can accurately identify modes in widely 
differing flight conditions without modification of 
algorithm parameters. 

(3) Certain pole-zero characteristics of the system are 
poorly identified and require interaction by the 
analyst to modify controlling parameters of the 
algorithm, such as assumed model order or frequency 
range. 

(4) In general, the number of modes assumed in the model 
should be equal to or exceed the number of modes 
actually present, for greatest accuracy. 

(5) The general "rules" described in Chapter III for 
selecting algorithm parameters give generally the 
best results with XV-15 data. 
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CHAPTER V 

APPLICATION TO CANTILEVER WING WIND TUNNEL DATA 

5.1 INTRODUCTION 

Tests were conducted by the government on a scaled canti- 
lever wing model in the M.I.T. wind tunnel for the purpose 
of measuring dynamic stability. Runs were made at four tunnel 
speeds and four inputs (8p, flap deflection; B^ s , longitudinal 
cyclic; 0 Q , collective and Gp c , lateral cyclic) were applied 
successively (see Table 5.1). At each speed and for each 
applied input, the output measurements are chordwise bending, 
q^, vertical bending, and the first two rotor flapping 

modes, (3p c and P^ s * In addition, runs were made without 
any applied input, where the model was excited by random 
tunnel turbulence. The sampling interval was 0.02 sec. The 
data was Fourier transformed and averaged on the basis of 
about 20 runs. The averaged autospectrum and cross spectrum 
of the inputs and outputs were supplied by the government for 
evaluation of the algorithm developed in Chapter III. A typical 
transfer function obtained by dividing the autospectrum of 
the input into the cross spectrum of the input and output 
are shown in Figure 5.1. In spite of the averaging, the 
transfer functions are quite noisy. The following sections 
describe the results. 

5.2 SINGLE INPUT, SINGLE OUTPUT 

The collective pitch input, 0 O , excites primarily the 
wing chordwise motion, and the flaperon input, 

excites primarily the wing vertical motion, q^, at all 
speeds. Therefore, these input-output pairs are used to eval- 
uate the algorithm. 
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Table 5.1 


Summary of Wind Tunnel Test Points 
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Several runs are made using input gp and output q^ 
to study the effect of various algorithm parameters on estima- 
tion accuracy. The highest frequency used for this input 
output pair is 10 Hz, and the model order is set to two (one 
DOF) . The minimum acceptable value of the delay parameter 
t is therefore (2x2)/(10x2) = 0.2 sec. With a x of 0.3 
and 7 < o) < 10, a two-pole, one-zero model gave the results 
shown in Table 5.2. Except when going from a nondimens ional 
speed of 0.47 to 0.58, the damping ratio increases with speed. 
If the same model is used except without any zeros, the esti- 
mates obtained are again shown in Table 5.2. The estimates 
when the delay parameter, t , is set equal to zero with and 


Table 5.2 

Estimate of the q^ Mode from Input 8p and Output q^ 
(Two Pole Model Assumed) 




V 

NATURAL 

DAMPING 

ZERO 

. 

GAIN 



nR 

FREQUENCY 

RATIO 

a. 

t = 0.3 

0.34 

8.29 

0.0167 

13.1 

-0.123 


7 < os < 10 

0.47 

8.23 

0.0212 

11.0 

-0.269 


0.58 

8.20 

0.0198 

9.54 

-0.446 


ONE ZERO 

0.69 

8.16 

0.0250 

12.2 

-0.579 

b. 

x = 0.3 

0.34 

8.29 

0.0157 



1.61 


7 < 03 < 10 

0.47 

8.23 

0.0192 

— 

2.90 


0.58 

8.21 

0.0182 

— 

4.21 


NO ZEROS 

0.69 

8.16 

0.0226 

— 

7.03 

c. 

7 < 03 < 10 

0.34 

8.30 

0.0108 

13.2 

-0.121 


x = 0.0 

0.47 

8.24 

0.0151 

11.0 

-0.268 


0.58 

8.21 

0.0148 

9.67 

-0.440 


ONE ZERO 

0.69 

8.16 

0.0195 

12.6 

-0.561 

d . 

7 < 03 < 10 

0.34 

8.27 

0.00656 

. 

1.53 


x = 0.0 

0.47 

8.19 

0.00723 

— 

2.78 


0.58 

8.16 

0.00642 

— 

4.08 


NO ZEROS 

0.69 

8.11 

0.0103 

— 

6.78 
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without a zero in the transfer function are also given in 
Table 5.2. A comparison of these results shows that the esti- 
mates of natural frequencies , zero locations and gains do not 
depend upon the value of z or the absence of zero from the 
transfer function. The estimate of the damping ratio , however, 
depends strongly on the correct choice of t and the number 
of zeros in the transfer function (see Fig. 5.2). 

The response of the output q£ to input B q is used to 
estimate the mode. A two-pole, one- zero model is used and 

t is set to 0.3. Two cases are tried with 5 < co < 15 and 
11 < a) < 15. The results are shown in Table 5.3 and Figure 5.3 
In all but one case, the algorithm indicates the presence of 
a zero in the transfer function. 



Figure 5.2 Estimate of the Damping Ratio of the q^ Mode 
Obtained from the Output q^ and Input S ^ 
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Table 5.3 

Estimate of the q£ Mode from Input 0 q and Output q 7 



V 

NATURAL 

DAMPING 

ZERO 

GAIN 


f2R 

FREQUENCY 

RATIO 

a. TWO POLE 

0.34 

12.9 

0.0396 

-137 

0.201 

ONE ZERO 

0.47 

12.5 

0.0327 

-18.3 

1.61 

7 < co < 15 

0.58 

12.4 

0.0250 

-17.5 

1.90 

t = 0.3 

0.69 

12.4 

0.0238 

-9.84 

1.09 

b. TWO POLE 

0.34 

12.9 

0.0347 

148 

-0.154 

ONE ZERO 

0.47 

12.6 

0.0308 

-45.6 

0.594 

11 < co < 15 

0.58 

12.5 

0.0269 

-43.9 

0.319 

t = 0.3 

0.69 

12.4 

0.0257 

-23.3 

0.555 
















5.3 SINGLE OUTPUT, MULTI -INPUT 


The algorithm is applied to measurements of q^ resulting 

from inputs and e ]_ 3 - Figure 5.4 shows the estimate of 

the damping ratio of the q-, mode based on the selection of 

1 V 

different model orders, for ^ = 0.47. The value of the 
time delay parameter is set to the minimum value corresponding 
to the model order selected in the estimation. It is clear 
from Figure 5.4 that the estimate of the damping is relatively 
independent of the model order. When the degrees of freedom 
are set to two or three, all but a pair of complex poles cancel 
out complex zero pairs. 


The effect of t on the estimate of the damping ratio is 
shown in Fig. 5 . 5 . Note that the estimate is essentially 
constant for t > 0.6 (the minimum acceptable value). Figure 
5.6 shows the estimate of q^ damping ratio at various speeds 
and compares the results with single input, single output cases 
presented in the previous section. 



Figure 5„4 Effect of Model Order on the Estimate 

of q^ Damping 
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ESTIMATED DAMPING RATIO OF THE 









5.4 MULT I -INPUT, MULTI -OUTPUT 


The responses of the measurements q.. , p and P^ 5 to 
inputs 6p and Q-^ c are used to estimate the modes of the 
cantilever wing. A six-pole, five-zero model is used and 
i is set at 0.6. The estimate of the damping ratio of the 
mode is shown in Figure 5. 7 for various speeds. These 
results agree quite closely with the case when only output 
q-^ is used. This indicates that P^ c and p^ s are not use- 
ful in estimating the mode. They may, however, be able 

to estimate other rotor modes. 



Figure 5.7 Estimate of q-^ Damping from Multi -input. 

Multi-output Data 
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5 . 5 SUMMARY 


The estimates of the various modes, obtained from the 
wind tunnel data through the application of the algorithm 
developed in Chapter III, are quite reasonable. The basic 
trends of actual frequencies and damping ratios agree with 
analytical results and results based on other techniques. 
Further processing of more wind tunnel data is required to 
establish the best procedure for using the algorithm. 
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CHAPTER VI 
CONCLUSIONS 


6.1 ALGORITHM DEVELOPMENT 


This study develops an on-line algorithm based on advanced 
system identification technology for real time measurement of 
dynamic stability in the Ames 40- by 80-foot Wind Tunnel. 

The algorithm is a four-step procedure requiring: (1) Fourier 

transforming the time history data, (2) setting up a set of 
linear equations, (3) solving these equations, and (4) computing 
roots of characteristic polynomials to find poles and zeros. 

The entire procedure is computationally efficient and appears 
to meet the requirements of the wind tunnel test. In addition, 
the algorithm can handle multi- input, multi-output data in a 
statistically optimal manner and minimizes the effects of pro- 
cess and measurement noise sources. Application of the algo- 
rithm to ground resonance simulation data indicated the follow- 
ing important features of the algorithm: 

(1) The algorithm provides good estimates of natural 
frequencies of all modes under every condition used 
in the simulation. The damping ratio estimates are 
excellent up to 20% random noise in the output, when 
a theoretically derived value of the time delay para- 
meter x is used in the algorithm. The estimates of 
the damping ratio degrade with increasing output 
noise and with nonoptimal t . 

(2) The accuracy of the estimates improves when more out- 
puts are used in identification. In fact, with the 
measurement of the lateral deflection alone, the rotor 
modes cannot be estimated. 

(3) The sampling interval and the data length are import- 
ant variables in determining estimation accuracy. 

The data record should be long enough to provide good 
frequency resolution in the region of the modes of 
interest. 
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(4) Discrete frequency noise (e.g., 1/rev noise) could 
degrade estimation accuracy significantly if it is 
not properly accounted for. However, the algorithm 
can be modified such that discrete frequency noise 
does not affect parameter estimates significantly. 

The application of the algorithm to the cantilever wing 
model indicated, in addition to the above, that: 

(1) If a high order model is specified in the estimation 
algorithm, several estimated poles are very close 
to estimated zeros. If pole-zero cancellation is 
performed, the correct number of poles and zeros 
are identified. 

(2.) The computation time increases with the model order 
specified in identification, the number of outputs 
and the number of inputs. The dependence on the 
model order is the strongest. 

(3) Scaling of inputs and outputs may be required for 
improved numerical conditioning. 


6.2 XV- 1 5 MODE IDENTIFICATION 


Applying the new algorithm to 13 DOF simulated XV- 15 
time history data leads to the following conclusions: 

(1) The algorithm achieves excellent results identifying 
the frequencies of significant aeroelastic XV-15 
modes, and generally very good results in the iden- 
tification of damping ratios. Cases of reduced 
accuracy in the estimated damping ratio are generally 
those in which multiple roots--poles and zeros-- 
occur at the same damped natural frequency in the 
basic system, a situation which changes with flight 
condition. 

(2) Good mode identifications are achieved over a wide 
range of flight conditions with no changes in the 
algorithm. Improvements can, in some cases, be 
obtained through specific changes made by the analyst 
in certain model parameters, such as assumed model 
order and frequency range, and in some cases, scaling. 
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(3) In general, assuming a model with a greater number 

of modes than actually anticipated gives much- improved 
accuracy in the identification of the most important 
modes, with "excess" roots being placed by the algo- 
rithm either well outside the frequency- damping 
range of interest or in self-cancelling pole- zero 
pairs. Attempting to use too low a model order 
generally gives poor damping predictions. 

(4) The benefits obtained from using multiple measure- 
ments appear to depend on the degree to which the 
measurements are related. With essentially independ- 
ent measurements, no improvement in identification 
was found using two measurements instead of one. 

C5) The basic "rules" on algorithm parameter specifica- 
tion developed on smaller models are valid for high- 
order models. Particularly careful attention must be 
paid to the effects of scaling. 


6.3 APPLICATION TO WIND TUNNEL DATA 


The algorithm of this study is successfully applied to 
data from the wind tunnel test of a cantilever wing model 
(Chapter V ) . The identification is performed off-line through 
on-line conditions are simulated as much as possible. The 
following conclusions may be drawn from this limited processing: 

(1) Zeros play an important role in accurate estimation 
of damping ratios. Neglecting a neighboring zero 
could significantly alter the identified value of the 
damping ratio. In addition, a proper selection of 
the time delay parameter r is important. 

(2) Estimates of damping ratios, from different sets of 
inputs and outputs, are close to each other. This 
verifies the results of the algorithm. 
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6.4 RECOMMENDATIONS 


Based on the studies completed thus far, the following 
recommendations are made regarding the direction of future 
development work: 

(1) Additional independence from judgments on the part 
of the analyst should be sought; e.g., selection of 
model order, scaling, etc. 

(2) Improved resolution of close pole- zero groupings 
should be sought. 

(3) Direct measurements of rotor degrees-of-freedom are 
important in many cases. Methods to extract these 
time histories from single blade measurements is 
required. 

(4) The algorithm must be further tested to derive opera 
tional procedures for effective use of the algorithm 

(5) The approach should be thoroughly tested on real 
wind tunnel data before it can be finally accepted 
as a method for flutter analysis. 
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APPENDIX A 

STATE, GENERALIZED STATE, AND ARMA MODELS 


This appendix discusses briefly the relationships between 
the generalized state and state variable formulation of the 
equations of motion, and shows that the state variable formula- 
tion is equivalent to the autoregressive moving average (ARMA) 
formulation used in the parameter identification method of 
this study. 

The generalized state equations are written in the form, 

A 0 x + A, x + A x = Bu 
2 lo 

where A„ , A-,, A n and B are coefficient matrices and x 

is a vector whose n elements are the degrees of freedom of 
the physical system. This set of equations may be transformed 
uniquely to state variable form by solving for x and defining 



hence , 


— 1 

o 

M 


o 

— 1 

_l 

X + 


r 

-A^A ! -A~ l A, 

2 0 i Z 1 -j 



(2n x 2n) (2nxq) 


= Fx + Gu 

The inverse transformation is, however, not unique. From the 
state variable form we obtain. 
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Cx + CAo^A, x + CAl^A x = CA^Bu, 

2 1 2 o 2 * 

where C is any nonsingular matrix. If it is known that 
C=A 2 , the original generalized state equation is obtained. 

The ARMA form describes the output, y, where 
y = Hx + Du 

Taking successive time derivatives of this expression yields 
y = Hx + Du = HFx + HGu + Du 
y = HFx + HGu + Dii 

= HF 2 x + HFGu + HGu + Du 

y(n) = Hf n^ + HF n-l Gu + H F n ' 2 Gu + ... + Du 1 ^ 

The characteristic polynomial of F is 

N (A) = A n + a n . 1 ^ n ' 1 + ... + a Q = 0 , 

hence, by the Cayley-Hamilton theorem, we may write 

F n = -a .F 11 " 1 - a ? F n " 2 - ... - a.F - a 
n-1 n-2 1 o 

Therefore , 

= - I a -HF 1 ^ + HF n ' 1 Gu + HF n ' 2 Gu +. . .+ Du^ 
i=0 1 

= Y a . (y (l) -HF 1_1 Gu-. . . -Du (l) ) +HF n " 1 Gu+ . . .+Du Cn} 
i=0 1 
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or 


y tn) + S = S 1 a - (HF 1_1 Gu + HF 1_2 Gu + . . .+ Du (l) ) 

i=0 1 i=0 1 


+ HF n-1 Gu +. . .+ Du^ 


= B + B -jU^ n_1 ^ + ... + B u. 

n n-i o' 


which is the ARMA model upon which the new algorithm was based. 
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APPENDIX B 

ON-LINE PARAMETER ESTIMATION USING THE 
INSTRUMENTAL VARIABLES APPROACH 


This appendix describes an algorithm for estimating the 
parameters of the autoregressive moving average model (Eq. 

(3.1) ) using the instrumental variables approach. Equation 

(3.1) may be written as 

y^ (t) + I a. I y( n_1 ) (t) = 
i = l 1 

Z (t) + (t)) (B.l) 

i=0 1 1 

Define 

T T 

U (t) = (u (n) (t), u Cn_1) (t), ... u T (t)} 


Y t (t) 

- {y> 


... YiCt)} 


a T - 

' ta l- a 2 

» V 



Cd 

H* 

II 

CCB 0 1 T . 

1 

(B ) T , ... (B ) T ) 
i i n i 


(B . 2) 

where y.(t) is the ith element of vector 

y(t) and 

(B,) is 
J i 

the ith column of B.. 

1 

Let 



■Y x (t) 

U(t) 0 ... 0 ' 



Q (t) = 

Y 2 (t) 

0 U(t) ... 0 


(B . 3) 


-V^ 

0 0 ... U(t)_ 
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Then Eq. (B.l) may be written as 


y^(t) = Q(t) 


B- 


\v 


+ Z (t) 

i=0 


CB.4) 


Premultiply both sides of the above equation by P(t,x) given by 


PCt,t ) = 


Y^Ct-x) Y 2 T (t-x) .. 


U T Ct) 


u T ct) 


0 

0 


V ct ’ T) 

. . 0 


. . 0 


u T Ct) 


CB.5) 


and take the expected value. For t>0, the expected value of' 


n 


P(t 5 x) 2 

i=0 1 


is zero. Therefore, 


E(PCt,x)y^^ Ct))= ECP(t,x)QCt)) 


'»V 

Bn 


Vi I 


(B.6) 


where E is the expected value operator. .Let 
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(B . 7) 


R ab C T ) ■ E(u(t+x)y T (t)) 
R T Ct) - tr(R aa (t)) 


and use the following identities 


(i + j ) 


(a) E(a [i) CtH-T)b ( ^ } (t)) = (-l) j R ab (T), x>0 


(B . 8) 


(b) R ab (x) = R ba (-t) 


The indicial superscript on R ab (x) implies differentiation with 
respect to t. Then 


n 


E(P(t,x)v y (t) A 


(B . 9) 


b = 
o 


(-D 


n-1 n (2n-l) 

iV- rj-» 

y y 


(x) 


C-l) n " 2 R^ n ' 2 ) (x) 

y y 


(-1)° R C ^ (t) 

y y 


n elements 


100 



c _ i ) n R^n) (0) 


R C2n-i) ^ 
uy. 1 ; 


1=1,2, . . .p 

Cn+1) elements 


R^CO) 


E(P(t,r)Q(t)) A 


A 11 

1 A 1 

I A 12 
l 

A 21 

1 A 

1 a 22 _ 

‘ v— ' 

■—v— ' ■ 

n 

Cn+1) pq 


} (n+'l) pq 


(-l) 11 " 3 ^ 211 ' 2 ) (r) C-l) n_1 R 

y y 


n ~ 1 R^ 11 ' Cx) 

y y 


(-l) n ' Z R (2n ' 3) (t) (-l) n ' 2 R C2n ' 4:) 

y y y y 


C-1)°R^'^(t) 

y y 


C-D°r^' 2) Ct) 

y y 


C-l) n " 1 R C ^" 1) Cx) 


Cx)-... C-l) n ’ 2 R C T" 2) Cx) 


... (- 1) U R T Cx) 


A 12 ^ A 12^ ’ ^ A 12^ . **• ( ^ A 12- ) * 
1 2 T3 


( A iPi - Fc-i ) 0 ; 1 r ^"' 15 ct) 


C-D° Ct) 


c-d "" 1 r^; 15 ct)‘ 


c- 1) ° R T Cx) 

UV . ^ J 


CB.14) 
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v 2 1 


( A 2i) 


(A 21^ 


( a 2 i) 


C B . 1 5 ) 


^ A 21^ . 

1 


r c-i) n ' 1 R^ n ' 1 ' ) co) 


n- 1 


Cn-1) 


L c-ir x R uy< co) 


. (- 1 ) ° R ^ n ' 1 ( 0 ) 

^ uy . v y 


(-D°R U)r .(0) 


(B . 16 ) 


22 


A 2 2 0 . . . 0 
0 


L 0 


l 2 2 . 

0 ' A 


0 

22 


(B.17) 


a 22 - rC-l) n R, C ,) n) C0) 


uu 




uu 


(-1) °R ( ^ (0) 
^ ' uu ^ } 


(-D u R uu C0) 


CB. 18) 


In the light of Eqs. (B.9) - (B.18), Eq. (B.6) may be written as 


b o = A ll a + A ^ A 1 2 ^ . 6 i 
i=l i 


(B . 19) 


b i = ^ A 21^. a + A 22 S i > i = l » 2 , . . . p 


(B. 20) 
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Equation (B. 20) may be solved for 6^ 


e i ~ A 2 2 ^ b i ' ( A Zl) . 

i 


(B. 21) 


Therefore, the equation for determining a is 

b o ‘ ^ A 12^ , A 2Z b i = ^ A ll " . z ^ A 12^ , A 2Z ^ A 21^ . }a 

x=l i 1=1 l i 


(B. 22) 


This is a set of n linear equations in n unknowns a. If the 
zeroes of the transfer function are desired they may be obtained 
by substituting a in Eq. (B.22). 

Estimation of Differentials of Correlation Functions: 


Matrices A's and vectors B's in Eq. (B.2 2) are functions of 
differentials and auto- and cross- correlations of input u and 
output y. These differentials must be estimated from data, which 
would typically be discrete input and output time histories (or 
equivalently the corresponding Fourier transforms) . There are 
two approaches which may be used to estimate differentials of 
the correlation functions. Let the inputs ’and outputs be u(i) 
and y(i) , i=l,2,...N, respectively. 


In the direct approach the auto- and cross -correlations of 
the inputs and the outputs are computed. 


yy 


uy 


uu 


CrA) = 

1 N-r 
K * 


i=l 


, N-r 

CrA) = 

I V 
M L 

lN i=0 

(rA) = 

, N-r 
1 >0 


.Tr: 


(B. 23) 
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where A is the sampling interval. r is chosen to be a multiple 
of A. The derivatives of the correlation functions are determined 
by numerical differentiation, e.g. 


RyyhO * ^ (R yy C^A) ' V CT ' 4) 

R^CO = (R yy (x»A) - R yy Ct-A) - 2R yy (x)) (B.24) 

etc . 


In the second approach, the differentials of correlation 
functions are estimated from the Fourier transforms of the inputs 
and the outputs given by 

Y (k) = S yCiOe’J 2 ^ 1 - 1 )^ 
i = l 

CB .25 

U(k) = Z uCi)e'^ 2lTCi - 1)k/N 
i=l 

k = 0,1,2, . .‘.N-l 


The autospectrum of y and the cross-spectrum of y and u are 
given by (for N even) 


S yy OO - I Y(10Y*0 c) 
S yu OO - | *0011*00 


The autocorrelation of y is given by 


R yy M 


N/2 - 1 j 2'TTxk 

L S yy (k)e 


(B. 26) 



1 


(B . 2 7) 
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Therefore , 




N/2 - 1 

y s (k) 

L, yy 




j 2 rk \ n ' j 2rrk (B.28) 


NA / NA 


Similarly , 


R^ co 


N/2 - 1 


v 

k=- l* l 


uy 


00 ( ; 


j 2 tt k \ n . j 2rtk 
NA / e NA 


(B. 29) 


and 


(x) 
uu ^ 


N/2 - 1 

y s rk) 

Lj uu ^ 

R=-| + l 


/i 2 7rrk \ 
j 2rk \ n \ NA / 
NA / e 


(B. 30) 


Equations (B.28) - (B.30) are very flexible. Since the system 
is linear, S (k) , S (k) , and S uu (k) may be passed through 
the same filter and then used in Eqs . (B.28) - (B.30). This 

is equivalent to passing the input and output data through the same 
filter before identification. It is useful when there is discrete 
frequency noise or a certain frequency range does not contain much 
information about the system modes. A band pass filter can then 
be applied to reject these frequencies. 
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APPENDIX C 


XV-15 MODEL AND DYNAMIC CHARACTERISTICS 


This appendix presents additional information on the XV-15 
model and on its predicted aeroelastic characteristics at the 
six flight conditions studied in this report (Chapter IV) . 

These data comprise the analytical solutions to which the 
results of algorithm application were compared. The contents 
of this appendix are: 

Table C.l: Analytical model degrees of freedom and modes. 

Table C.2: State, control, gust, and observation vectors. 

Table C.3: Wingtip deflection constants. 

Tables C.4 to C.9: Roots of symmetric wingtip deflection/ 

control input transfer functions. 

Tables C.10 to C.15: Roots of antisymmetric wingtip 

deflection/ control input transfer functions 




Y/Oj 1 
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Table C.l 

Analytic Model Degrees of Freedom and Modes 


DEGREES OF FREEDOM 

MODES 

SYMMETRIC 

ANTISYMMETRIC 

SYMMETRIC 

ANTISYMMETRIC 

B o 


S o 

B o 

6 lc 

B lc 

6 +l 

e n 

B ls 

e is 

e -l 

B -i 

he 

^lc 

? + l 

? +l 

? ls 

c ls 

S -1 

? -l 





a 

a 

a 

a 

<1 

q l 

q i 

q i 

q 2 

q 2 

q 2 

q 2 

P 

P 

P 

P 

BV 

BV 

BY 

BY 

S'L 

SS 

SL 

SS 

BL 

BS 

BL 

BS 
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T ab 1 e C . 2 

State, Control, Gust, and Observation Vectors 


Symmetric " ^ic c is 6 o 8 ic 6 is BL q i SL q 2 p BV 

X ANTI- = feic ? is B o B lc B ls BY BS q 2 q l SS p 

SYMMETRIC 

X T = [X T i X T ] (SYMMETRIC AND ANTISYMMETRIC) 

y T = [X TIP Z TIP ] 

= f 6 o 6 lc 9 ls 6 F^ 

g T = [u g V g w g ] 

T T 1 T 

u' = [u l j g 1 ] 

NOTE: Vector elements are defined in "Nomenclature", pp. vi-x. 
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Table C . 4 

Symmetric Transfer Function Roots 

Flight Condition #1: Hover 

a p = 90° 

V - 0 

n = 565 RPM 


MODE 

POLES 

ZEROS 

Z TIP /9 o 


REAL 

IMAGINARY 

REAL 

IMAGINARY 

s 

-0.023 

0 

0 

0 


-0.071 

0 

-0.049 

0 

BL 

-0.0021 

+0.202 

-0.0021 

+0.202 


-0.0087 

+0.286 

-0.0070 

+0.348 

SL 

-0.0067 

+0.348 

-0.264 

+0. 107 

2-1 

-0.264 

+0.108 

-0.089 

+0.465 

5-1 

-0.089 

+0.467 

-0.0077 

+0.573 

^2 

-0.0082 

+0.574 

-0.0091 

+0.751 

P 

-0.0100 

+0.738 

-0.0080 

+0.817 

BV 

-0.0083 

+0.820 

-0.0146 

+0.933 

a 

-0.0132 

+0.952 

-0.316 

+1.91 

6 , 

-0.316 

+1.91 

-0.0180 

+2.53 






B 

-0.225 

+2.21 

-0.426 

+16.6 

0 





5+1 

-0.0155 

+2.54 




111 















Table C.5 

Symmetric Transfer Function Roots 


Flight Condition #2 
a = 82.3° 

r 

V = 80 kts 
8 = 565 RPM 


MODE 

POLES 

ZEROS 

W 9 a 

z tip /6 f 


REAL 

IMAGINARY 

REAL 

IMAGINARY 

REAL 

IMAGINARY 


-0.027 

±0.022 

-0.100 

0 

-0.0282 

±0.020 

BL 

-0.0022 

±0.202 

+0.073 

0 

-0.0019 

±0.202 

ii 

-0.0129 

±0.291 

-0.192 

±0.093 

0.0002 

±0.325 

5L 

-0.0080 

±0.353 

-0.0026 

±0.213 

-0.262 

±0.116 

6 -i 

-0.2659 

±0.106 

-0.0092 

±0.319 

-0.0725 

±0.542 

? -i 

-0.0751 

±0.549 

-0.0746 

±0.538 

0.0472 

±0.617 

q 2 

-0.0102 

±0.575 

-0.0146 

±0.631 

-0.0764 

±0.650 

P 

-0.0091 

±0.737 

-0.0095 

±0.700 

-0.0861 

±0.807 

BV 

-0.0089 

±0.820 

+0.0373 

±0.844 

0.0621 

±0.819 

a 

-0.0142 

±0.956 

-0.0809 

±0.857 

-0.311 

±1.909 


-0.307 

± 1.92 

-0.363 

±1.81 

-0.262 

±2.235 

S o 

-0.260 

±2.26 

-0.070 

±2.57 

-0.0085 

±2.653 


-0.008 

±2.63 

2.77 

0 

— 

— 




-3.81 

0 
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Table C .6 

Symmetric Transfer Function Roots 

Flight Condition #3 
a p =60° 

V = 100 kts 
fl = 565 RPM 


MODE 

POLES 

ZEROS 

=w e „ 

W S F 


REAL 

IMAGINARY 

REAL 

IMAGINARY 

REAL 

IMAGINARY 

*3 

-a.019 

±0.028 

0.0116 

0 

-0.020 

±0.029 

BL 

-0.0023 

±0.203 

-0.0162 

0 

-0.0017 

±0.201 

‘h 

-0.0097 

±0.292 

-0.0018 

±0.231 

+0.0064 

±0.335 

SL 

-0.0098 

±0.397 

-0.0053 

±0.284 

-0.310 

±0.135 


-0.310 

±0.126 - 

-0.312 

±0.121 

-0.0118 

±0.464 

c -i 

-0.062 

±0.496 

-0.059 

±0.488 

-0.072 

±0.498 

^2 

-0.0106 

±0.573 

-0.0193 

±0.556 

-0.0116 

±0.723 

P 

-0.0087 

±0.743 

-0.0088 

±0.742 

-0.0092 

±0.740 

BV 

-0.0094 

±0.814 

-0.0115 

±0.843 

-0.0135 

±0.947 

a 

-0.0134 

±0.938 

-0.Q20 

±0.907 

-0.296 

±1.89 

Ai 

-0.298 

±1.89 

-0.272 

-1 .87 

-0.276 

±2.25 

6 

0 

-0.274 

±2.26 

-0.319 

±2.58 

-0.0164 

±2.61 

C +1 

-0.0128 

±2.58 

4.55 

0 

— 

— 




-4.77 

0 
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Table C.7 

Symmetric Transfer Function Roots 


Flight Condition #4 

a p = 30 ° 

V = 120 kts 
57 = 565 RPM 


MODE 

POLES 

ZEROS 

X TIP /0 o 

z tip /5 f 


REAL 

IMAGINARY 

REAL 

IMAGINARY 

REAL 

IMAGINARY 


-0.0320 

±0.031 

+0.0073 

0 

-0.0336 

• ±0.032 

Bl 

-0.0026 

±0.202 

-0.0091 

0 

-0.0022 

±0.203 

Ol 

-0.0056 

±0.296 

-0.0021 

±0.239 

-0.0535 

±0.35 

«-l 

-0.0645 

±0.353 

-0.0044 

±0.302 

-0.0034 

±0.426 

St 

-0.0101 

±0.435 

-0.0638 

±0.345 

-0.325 

±0.164 

B -1 

-0.3157 

±0,138 

-0.319 

±0.131 

-0.0375 

±0.652 


-0.0095 

±0.566 

-0.0096 

±0.568 

-0.0387 

±0.677 

P 

-0.0079 

±0.680 

-0.0110 

±0.699 

0.364 

±1.15 

BV 

-0.0095 

±0.822 

-0.0112 

±0.837 

-0.683 

±1.09 

a 

-0.0141 

±1.004 

-0.0186 

±0.988 

-0.200 

±1.81 

B +1 

-0.320 

±1.88 

-0.2847 

±1.87 

-0.262 

±2.16 

B o 

-0.253 

±2.21 

-0.0194 

±2.43 

1.11 

±3.07 

*+1 

-0.0216 

±2.44 

-0.0809 

±8.30 

— 

— 
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Table C.8 

Symmetric Transfer Function Roots 
Flight Condition #5 

a p " °° 

V = 160 kts 
Q = 458 RPM 


MODE 

POLES 

ZEROS 

X TIP /e o 

wv 


REAL 

IMAGINARY 

REAL 

IMAGINARY 

REAL 

IMAGINARY 

> 

s 

-0.0646 

±0.016 

0 

0 

-0.0646 

±0.016 

■ 

BL 

-0.0031 

±0.250 

-0.0016 

0 

-0.0030 

±0.249 

q i 

-0.0096 

±0.372 ‘ 

-0.0027 

±0.296 

-0.309 

±0.136 

c -i 

-0.0933 | 

±0.405 

-0.0096 

±0.371 

-0.092 

±0.394 


-0.293 i 

±0.142 

-0.094 

±0.406 

+0.0079 

±0.587 

SL 

-0.0117 

±0.560 

-0.292 

±0.141 

-0.0176 

±0.627 

P 

-0.0205 

±0.809 

-0.165 

±0.776 

-0.0176 

±0.840 

^2 

-0.0X11 

±0:838 

-0.0103 

±0.816 

-0.0021 

±1.25 

BV 

-0.0124 

±1.03 

-0.0119 

±1.02 

-0.0035 

±1.51 


-0.316 

±1.84 

-0.306 

±1.82 

-0.290 

±1.89 

a 

-0.013 

±2.04 

-0.020 

±1.96 

-0.094 

±2.48 

6 o 

-0.066 

±2.55 

-0.0607 

±2.55 

-0.209 

±2.57 


-0.201 

±2.57 

-0.104 

±5.23 

— 

— 
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Table C.9 

Symmetric Transfer Function Roots 


Flight Condition #6 

«p " °° 

V = 200 kts 
R = 458 RPM 


MODE 

POLES 

ZEROS 

X TIP /6 o 



REAL 

IMAGINARY 

. REAL 

IMAGINARY 

REAL 

IMAGINARY 


-0.041 

0 

0 

0 

-0.041 

0 


-0.152 

0 

-0.0032 

0 

-0.152 

0 

BL 

-0.0030 

±0.250 

-0.0026 

±0.296 

-0.0028 

±0.25 


-0.0897 

±0.314 

-0.090 

±0.315 

-0.0884 

±0.308 

**1 

-0.0111 

±0.372 

-0.0118 

±0.371 

-0.346 

±0.188 

6 -l 

-0.315 

±0.196 

-0.315 

±0.194 

+0.035 

±0.576 

SL 

-0.0112 

±0.558 

-0.0245 

±0.764 

-0.036 

±0.602 

P 

-0.0280 

±0.791 

-0.0076 

±0.809 

-0.020 

±0.839 

q 2 

-0.0102 

±0.835 

-0.0123 

±1.02 

+0.0038 

±1.25 

BV 

-0.0130 

±1.03 

-0.337 

±1.77 

-0.052 

±1.50 

6 ±1 

-0.351 

±1.79 

-0.0211 

±1.95 

-0.317 

±1.87 

a 

-0.0118 

±2.03 

-0.0768 

±2.47 

-0.104 

±2.39 

6o 

-0.0792 

±2.47 

-0.108 

±5.15 

-0.206 

±2.51 


-0.198 

±2.51 

— 

— 

— 

— 
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Table C.10 

Antisymmetric Transfer Function Roots 

Flight Condition #1 

a = 90° 

P 

V = 0 

fi = 565 RPM 



POLES 

ZEROS 

MODE 

Z TIP /e o 


REAL 

IMAGINARY 

REAL 

IMAGINARY 

BS 

-0.0032 

±0. 281 

-0.0012 

±0.211 

BY 

-0.0027 

±0.228 

-0.0018 

±0.240 

q 2 

-0.0062 

±0.330 

-0.0054 

±0.332 

q l 

-0.0114 

±0.391 

-0.269 

±0.113 


-0.0877 

±0.467 

-0.0964 

±0.462 

Kl 

-0.264 

±0.108 

-0.148 

±0.686 

h. 

-0.0196 

±0.577 

-0.0004 

±0.699 

p 

-0.0116 

±0.621 

+0.150 

±0.683 

ss 

-0.0088 

±0.707 

-0.0170 

±0.999 

a 

-0.0113 

±0.908 

. -0.312 

±1.93 

6 +l 

-0.318 

±1.914 

-0.075 

±2.54 


-0.211 

±2.364 

-2.54 

0 

q +l 

-0.0152 

i 

±2.53 

2.61 

0 
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Table C.ll 

Antisymmetric Transfer Function Roots 

Flight Condition #2 
a p = 82.3° 

V = 80 kts 
n = 565 RPM 


MODE 

POLES 

ZEROS 

X TIP /0 o 

z tip /6 f 


REAL 

IMAGINARY 

REAL 

IMAGINARY 

REAL 

IMAGINARY 

BS 

-0.0034 

±0.281 

-0.0069 

±0. 237 

-0.0031 

±0.229 

BY 

-0.0030 

±0.229 

-0.0030 

±0.281 

-0.0037 

±0.259 

q 2 

-0.0067 

±0.333 

-0.0108 

±0.312 

-0.0064 

±0.335 

q i 

-0.0140 

±0.397 

-0.238 

±0.088 

-0.272 

±0.109 

6 -i 

-0.2670 

±0.106 

-0.0081 

±0.538 

-0.0794 

±0.542 

s-i 

-0.0740 

±0.545 

-0.0784 

±0.544 

-0.0168 

±0.572 

*s 

-0.0158 

±0.578 

-0.0081 

±0.631 

-0.0265 

±0.690 

p 

-0.0099 

±0.630 

-0.0098 

±0.704 

+0.0135 

±0.711 

SS 

-0.0085 

±0.704 

-0.0051 

±0.939 

-0.0059 

±0.861 

a 

-0.0139 

±0.909 

-0.346 

±1.89 

-0.3100 

±1.92 


-0,307 

±1.91 

-0.0047 

±2.61 

-0.230 

±2.45 

6 o 

-0.239 

±2.45 

5.68 

0 

-0.0773 

±2.59 


-0.0101 

±2.62 

-6.19 

0 

-— 

— 
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Table C.12 

Antisymmetric Transfer ' Function Roots 
Flight Condition #3 

a p = 60 ° 

V = 100 kts 
a = 565 RPM 


MODE 

POLES 

ZEROS 

X TIP /e 0 

Z TIP /l5 F 


REAL 

IMAGINARY 

REAL 

IMAGINARY 

REAL 

IMAGINARY 

BS 

-0.0032 

±0.282 

-0.0032 

±0.280 

-0.0030 

±0.230 

BY 

-0.0031 

±0.231 

-0.0031 

±0.298 

-0.0028 

±0.270 

q 2 

-0.0053 

±0.352 

-0.0095 

±0.363 

-0.0066 

±0.361 ’ 


-0.0134 

±0.398 

-0.0616 

±0.491 

-0.0611 

±0.494 


-0.0654 

■ ±0.494 

-0.315 

±0.093 

-0.313 

±0.121 

h 

-0.0079 

±0.585 

-0.0032 

±0.628 

-0.0402 

. ±0.589 

6-1 

-0.309 

±0.126 

-0.0097 

±0.665 

0.0194 

±0.583 

SS 

-0.0110 

±0.663 

-0.0046 

±0.838 

-0.0103 

±0.658 

P 

-0.0093 

±0.702 

-0.261 

±1.47 

-0.0161 

±0.918 

a 

-0.0138 

±0.918 

0.798 

±1.34 

-0.300 

±1.89 

6.1 

-0.299 

±1.89 

-1.57 

±1.63 

-0.267 

±2.41 

6 o 

-0.262 

±2.43 

-0.0453 

±2.63 

0.0084 

±2.56 

«+l 

-0.0135 

±2.58 

— 

— 

— 

— 
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Table C.13 

Antisymmetric Transfer Function Roots 
Flight Condition #4 

a p = 30 ° 

V = 120 kts 
0! = 565 RPM 


MODE 

POLES 

ZEROS 

X TIP /e o 

z tip /5 f 


REAL 

IMAGINARY 

REAL 

IMAGINARY 

REAL 

IMAGINARY 

BY 

-0.0038 

±0.231 

-0.0024 

±0.266 

-0.0036 

±0.231 

BS 

-0.0031 

±0.282 

-0.0030 

±0.283 

-0.0020 

±0.272 

q 2 

-0.0073 

±0.375 

-0.0613 

±0.355 

-0.0627 

±0.352 

q l 

-0.0177 

±0.392 

-0.0083 

±0.378 

-0.0081 

±0.376 

«-l 

-0.0060 

±0.357 

-0.3111 

±0.129 

-0.316 

±0.132 

S -1 

-0.3148 

±0.137 

0.022 

±0.631 

-0.0142 

±0.620 

*s 

-0.0146 

±0.588 

-0.0288 

±0.648 

-0.0794 

±0.634 

SS 

-0.0140 

±0.653 

-0.0085 

±0.784 

0.053 

±0.629 

P 

-0.0086 

±0.759 

-0.0254 

±0.914 

-0.0228 

±0.993 

a 

-0.0171 

±0.944 

-0.293 

±1.96 

-0.322 

±1.88 

e + i 

-0.321 

±1.88 

-0.0201 

±2.44 

-0.272 

±2.30 

B o 

-0.267 

±2.32 

-11.65 

0 

-0.0057 

±2.44 

r , 

-0.0216 

±2.44 

16.1 

0 





’+1 








120 



















Table C. 14 

Antisymmetric Transfer Function Roots 
Flight Condition #5 

a p = 0 

V = 160 kts 
Q = 458 RPM 


MODE 

POLES 

ZEROS 

X TIP /e o 

Z TIP /e F 


REAL 

IMAGINARY 

REAL 

IMAGINARY 

REAL 

IMAGINARY 

q 2 

-0.0152 

±0.486 | 

0 

±0.184 

-0.0053 

±0.285 

BS 

-0.0041 

±0.348 

-0.0040 

±0.348 

-0.0041 

±0.339 

q i 

-0.0177 

±0.520 

-0.0938 

±0.401 

-0.0931 

±0.402 

BY 

-0.0044 

±0.284 

-0.295 

* ±0.139 

-0.0122 

±0.486 

5 -i 

-0.0929 

±0.402 

-0.0074 

±a.5io 

-0.293 

±0.138 

e -i 

-0.293 

±0. 140 

-0.0136 

±0.820 

-0.0576 

±0.729 

*5 

-0.0624 

±0.723 

-0.114 

±1.03 

-0.152 

±0.773 

SS 

-0.0183 

±0.820 

-0.0934 

±1.03 

+0.110 

±0.789 

P 

-0.0105 

±0.984 

-0.314 

±1.89 

-0.305 

±1.87 

6 +l 

-0.320 

±1.87 

-0.0781 

±2.44 

-0.0236 

±2.47 

^1 

-a. 0375 

±2.42 

-0.0559 

±2.80 

-0.237 

1 

±2.64 


-0.240 

±2.65 

-8.39 ; 

0 

-0.117 1 

| 

±3.36 

a 

i 

-0.0405 

±2.87 

8.26 

0 

! 

— 

— 













Table C . 15 

Antisymmetric Transfer Function Roots 


Flight Condition #6 

a p = 0 

V = 200 kts 
Q = 458 RPM 


MODE 

POLES 

ZEROS 

X TIP /O o 

z tip /6 f 


REAL 

IMAGINARY 

REAL 

IMAGINARY 

REAL 

IMAGINARY 

*■*2 

- 0.0142 

± 0.485 

0.0002 

± 0.196 

- 0.0057 

± 0.285 

q l 

- 0.0224 

± 0.517 

- 0.0906 

± 0.314 

- 0.0027 

± 0.341 

BS 

- 0.0042 

± 0.348 

- 0.0041 

± 0.349 

- O .0904 

± 0.312 

BY 

- 0.0044 

± 0.283 

- 0.317 

± 0.194 

- 0.0117 

± 0.484 

^-1 

- 0.0908 

± 0.314 

- 0.0072 

± 0.505 

- 0.316 

± 0.189 

6 -i 

- 0.317 

± 0.192 

- 0.0175 

± 0,813 

- 0.0799 

± 0.726 


- 0.0950 

± 0.716 

- 0.146 

± 1.03 

- 0.168 

± 0.759 

SS 

- 0,0234 

± 0.812 

0.116 

± 1.04 

0.102 

± 0.784 

P 

- 0. 0117 

± 0.980 

- 0.358 

± 1.87 

- 0.339 

± 1.84 

6 + l 

- 0.362 

± 1.84 

- 0.085 

± 2.34 

- 0.034 

± 2.38 

C +1 

- 0.0398 

± 2.33 

- 0.0527 

± 2.80 

- 0.232 

± 2.56 

6 o 

- 0.235 

± 2.56 

-8.22 

0 

- 0.118 

± 3.35 

a 

- 0.0423 

±2.86 

8.08 

0 

— 

— 
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